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Abstract 

In analyzing neutrino oscillation experiments it is often assumed that while 
new physics contributes to neutrino masses, neutrino interactions are given 
by the Standard Model. We develop a formalism to study new physics effects 
in neutrino interactions using oscillation experiments. We argue that the 
notion of branching ratio is not appropriate in this case. We show that a neu- 
trino appearance experiment with sensitivity to oscillation probability P^ p 
can detect new physics in neutrino interactions if its strength Gn satisfies 
(Gn/Gf) 2 ~ Pij P - Using our formalism we show how current experiments 
on neutrino oscillation give bounds on the new interactions in various new 
physics scenarios. 
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I. INTRODUCTION 



The goal of neutrino oscillation experiments is to probe those extensions of the Standard 
Model (SM) which predict non-vanishing neutrino masses. However, in the usual treatment 
of neutrino oscillation experiments, it is often assumed that neutrino interactions are de- 
scribed just by the SM 0. While we know that this is a good approximation, often physics 
beyond the SM induces also new neutrino interactions. If New Physics (NP) contributes 
significantly to neutrino interactions, the conclusions that we draw from the experimental 
data can be affected. For example, even for massless neutrinos, the NP can allow for weak 
eigenstate muon neutrino to produce an electron in the detector, and in this case, we may 
erroneously conclude that oscillations have occurred. 

To search for NP effects in massive fermion interactions (quarks and charged leptons) 
the experimentally measured branching ratios are used in a straightforward way. However, 
in the case of neutrinos, there are two important subtleties: 

• Neutrino masses are unknown, and their difference may be very small. In such a case, 
experiments cannot observe neutrinos as mass eigenstates, and the results are sensitive 
to the time evolution of the flavor eigenstates. 

• The neutrino flavor is identified by charged current interactions. Since NP may modify 
them, the identification of the neutrinos cannot be done in a model independent way. 

The results of neutrino oscillation experiments are sensitive to the following three ingredients: 
The production process, the time evolution and the detection process. It is impossible to 
separate the NP contributions to the neutrino production or detection process: a formalism 
that combines all the three ingredients is necessary. 

II. FORMALISM 

For simplicity, and without loss of generality, we assume two neutrino flavors, CP con- 
servation and that the neutrinos are highly relativistic. First, we define the bases we use. 
Since the mass basis is well defined, we express all neutrino states as superpositions of mass 
eigenstates. We always work in the mass basis for the charged leptons. Then, the weak 
interactions define the weak basis where the neutrino weak eigenstates are 577(2) partners 
of the charged leptons. We start with a specific example and generalize it later. We consider 
a muon neutrino beam produced by 7r — ► \w decay, and the subsequent detection of electron 
neutrinos through vn — > ep. The weak eigenstate is given by a superposition of mass 

eigenstates \v™) as 

l<) = E^IO, (i) 

a 

so that | | 2 oc \(i>™, fi\?i w \iT)\ 2 where 7i w is the weak interaction Hamiltonian. In the 
presence of NP there might be extra carriers of the charged current interaction besides the 
W boson. Therefore, the neutrino produced by tt — > \w may be different from \vv[). We 
define this neutrino as a source basis eigenstate given by a different superposition of 
mass eigenstates 
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K> = E^JO, (2) 

a 

so that |f/^| 2 oc \(v™,fi\H\ir)\ 2 , H = H w + H NP where H NP is the NP interaction Hamil- 
tonian. Similarly, we define the neutrino detected by vn — ep as a detector basis eigenstate 
given by another superposition of mass eigenstates 

k d > = E^JO, (3) 

a 

so that \Uf a \ 2 oc |(C,^|e,P>| 2 . 

In general, for any neutrino oscillation experiment it is useful to use the following bases: 

• The mass basis, where the neutrino mass matrix is diagonal. 

• The weak basis, ||^)|, where the leptonic couplings of the W are diagonal. 

• The source basis, {1^)}, where the interaction of the production process is diagonal. 

• The detector basis, ||^)|, where the interaction of the detection process is diagonal. 

When neutrino interactions are fully described by the SM, the last three definitions coincide, 
and this basis is usually called the interaction or the flavor basis. However, the main lesson 
from the above discussion is that in the presence of NP those three bases can be different. 

The source and the detector bases are related to the mass basis through the unitarity 
transformations 

k s > = E^IO, \4) = E^alO, (4) 

a a 

with £ — e, /i, t. The amplitude for finding a z/[ in the original v\ beam at time t is 

«WI)(t) = J2K\uZ e ~ tEatu X) = E e- iEat u s ea u? a , (5) 

where in the last step we have used the orthogonality of the mass eigenstates. The probability 
of finding a v% in the original v\ beam at time t is 



= E \ULU%U£U$f\ cos [{E a - Ep)t - arg (u^U^U d pn 

For two neutrino flavors and with CP conservation we have 

= / cos 9 ms - sin 6 ms \ = ( cos 9 md - sin 6 md \ 

V sin 9 ms cos 9 ms J ' V sin 9 md cos 9 md J ' 

with \9 ms \,\9 md \ < 7r/4. Define 



(6) 



(7) 
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x = ™ t , Am 2 = m\- m 2 , 9 sd = 9 md - 9 ms . (8) 

4£/ 



Using E a — Ep ~ (m 2 — m%)/2E, we get our main result: 

Pe/^z) = sin 2 9 sd + sin 20 md sin 29 ms sin 2 x. (9) 



Few points are in order: 

1. When neutrino interactions are described by the SM, 9 md = 9 ms = 9, and Eq.(||) 
reduces to the known result, P ef i(x) = sin 2 29 sin 2 x 

2. NP that affects the production and the detection processes in the same way cannot 
be detected in appearance experiments. In those cases the flavor eigenstates are the 
same for all processes, even if they may differ from the weak eigenstates. Then we 
have, 9 md = 9 ms and Eq.(^) reduces again to the standard form, so that we cannot 
distinguish this situation from the SM interaction case. 

3. Experiments are working with neutrino beams, not necessarily monoenergetic. P^ p 
is defined to be the total appearance probability of a specific experiment, where the 
dependence on the energy spectrum of the beam and on the baseline length L, is 
included. If neutrinos are produced by several decays of different initial states, then 

^7 = £ a ^7(*)' ( 10 ) 



where ai is the relative weight of the i'th decay mode in the neutrino beam, and P^iS) 
is the appearance probability had only the z'th decay mode been responsible for the 
neutrino production. 

4. In two limits, Am 2 ^> E/L (x — > oo) and Am 2 = (x — 0), P efl is x independent. 
When x is large, sin 2 x averages to |, then Eq.(||) gives 

P efl = sin 2 9 sd + - sin 29 ms sin 29 md . (11) 

For massless (or degenerate) neutrinos, x = and the appearance probability becomes 

P ef , = sm 2 9 sd . (12) 



We learn that a signal can be seen in appearance experiments even for massless neu- 
trinos. This is the case when 9 sd ^ 0, namely, when the interaction in the produc- 
tion process is different from the interaction of the detection process. This signal 
is constant in distance, and does not have an oscillation pattern. We conclude: a 

distance-independent signal is not enough to prove that neutrinos are massive. Only 
an oscillation pattern provides a proof. 
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Experimentally, we know that neutrino interactions are dominantly those of the SM. 
Therefore, while 9 ms and 8 m d may be large, their difference has to be small. It is therefore 
reasonable to work in the weak basis and treat the NP as a perturbation. In the two 
generation case we define two small angles, 9w s and Owd, that parameterize the deviation 
from the weak basis, and a third angle (not necessarily small) 9wm, that rotates from the 
weak to the mass basis. Using 9 a b + 6bc = 9 ac and 6 a b = —Oba, we get the relations between 
these angles and those defined in (|7p: 6 S( i = Owd — &Ws and 9 m d = 6wd — Ov/m- 

To find the rotation angles we use the previously mentioned example, but the final 
result is general. We consider a muon neutrino beam produced by n — > [iv decay, while 
the subsequent electron neutrinos are detected through vn —>■ ep. For SM interactions the 
produced neutrinos are vY[ , and the detected ones are . We are interested in NP that 
gives effective couplings of the form G s N udJiuY and G d N udev^ with \G S N \ ^ \G d N \. Then, the 
produced and detected neutrinos are superpositions of weak eigenstates, i/* ~ GpuJ^ -\-G s N vY 
and v d ~ G F v™ + G%v™ , and it follows that 



QS Qd 

Ows ~ 7^—, Owd ~ 7^—- (13) 



G1 " vv a (~i 



F 



It is useful to express the appearance probability in terms of the rotation angle from the 
weak to the mass basis, 9w m , and the NP strength, G S N and G d N . For i = 0we get from (|12|) 




- • (14) 



Experimentally we know that in the x — > oo limit all the relevant angles are small. Then 
we get from ( |TT|) 

P - W {^t) 2 + (S) 2 + Wwm { qk GV k ) + (15) 

From the above formulae we can obtain an order of magnitude estimate of the strength 
of the experimentally relevant NP. An experiment with sensitivity to oscillation probability 
P^ p can probe NP in neutrino interactions when its effective strength Gn ~ ma,x(G s N , G d N ) 
satisfies 



G p 



>p: x j>. (i6) 



III. BRANCHING RATIO 

Measurements of Branching Ratios (BR) are widely used in searching for NP effects. The 
meaning of BR is unambiguous when discussing quarks and charged leptons, for which a BR 
measures the transition rate between mass eigenstates. The main advantage of using BR 
is that only the production process is relevant to the calculation, and one need not worry 
how the decay products are detected. Therefore, measuring a BR has to be independent of 
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the experimental setup and of the theoretical model under study. Since experiments cannot 
detect neutrinos as mass eigenstates, for neutrinos these requirements are not satisfied: the 
time evolution of the flavor eigenstates, and the NP effects in the detection process cannot, in 
general, be separated from the analysis of the experimental results. Therefore, the extension 
of the BR notion to the neutrino case is problematic. 

We see three major disadvantages in using BR for neutrinos: 

1. Using BR calculations we can probe only part of the parameter space. In the two gener- 
ation case three parameters describe the results of the neutrino oscillation experiments 
(see Eq.(|)): Am 2 , 9 ms and 9 m d. The BR calculation is sensitive only to one parameter, 
the mixing angle that rotate from the source basis into the basis we are interested in. 
For example, the BR into final mass eigenstate is given by BR(7r — > /iz/j™) = sin 2 # ms . 

2. In order to compare BR calculations with experiments, the dependence of the experi- 
mental results on Am 2 and 8 m d has to be removed. However, this cannot be done in a 
model independent way since each kind of NP may contribute differently to neutrino 
masses and to the detection process. Therefore, experimental measurements of BR 
cannot be presented in a model independent way. 

3. Finally, there is a problem of definition. Are the theoretical calculations and the 
experimental bounds on rare decays as BR(7r — > fiu e ) @||, BR(K — > [iv e ) [4| and 
BR(/i — > ev e v^ ||§,f| related to the same quantities and therefore directly compara- 
ble? Calculations were done for neutrinos in the weak basis 0, and in the mass basis 

Experimental results are presented as bounds on the relative appearance proba- 
bility, Pej5 p {i) = P^ p /ai, where Oj is the relative weight of the relevant decay mode in 
the neutrino beam. In the x — > limit they correspond to bounds on the BR for neutri- 
nos in the detector basis, e.g. BR(7r — > fiuf). In the x — > oo limit, all the relevant an- 
gles are small, and v™ (v™) couples mainly to the electron (muon). Then, to first order 
in small angles, electrons are detected when is produced at the source, or when z/™ 
creates an electron at the detector. Therefore, the experimental results bound the sum 
of the rare BR in the production process and the ratio of the cross sections in the de- 
tecting process, e.g. BR exp (Tr — > fiu e ) = BR(n — > /i^™) + a{i>™n — > ep)/a(u" l n — > ep). 
We see that the calculations cannot be directly compared with the published experi- 
mental bounds. 

We conclude: The notion of BR can be used to probe only part of the parameter space. 
The BR cannot be measured in a model independent way, and comparisons of experimental 
and theoretical results have to be done always very carefully, paying attention to check that 
the definitions are consistent. Therefore, the notion of BR is not appropriate when searching 
for NP in neutrino oscillation experiments, and the formalism that we have developed here 
is preferable. 

IV. EXAMPLES 

We give now three examples of NP scenarios with non-standard neutrino interactions 
that can be probed by current and near future neutrino oscillation experiments. In each 
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case we first briefly present the model and the new neutrino interaction, then we discuss 
the actual experiment for probing the new interaction. Then, we find the neutrino sector 
independent experimental constraints on the strength on the new interaction, and show 
how results of neutrino oscillation experiments can put stronger bounds on it in part of the 
neutrino sector parameter space. For simplicity, we do not specify the NP responsible for 
neutrino masses. 

In the first example we consider the Minimal Supersymmetric SM (MSSM) without R- 
Parity |10|,|lTJ. We consider a very simple case where the MSSM superpotential is extended 
with only one extra term, X 12 3[L\, Lj^]E R , where L % L (E R ) are lepton doublet (singlet) su- 
per mult iplets. Via the exchange of the singlet charged scalar, f R = E R , such a term gives 
rise to the effective four fermion interaction (in the weak basis) fll] 



C = G N m lx P L e) (pJi X PlvY) + h.c. , (17) 

where Pl = (1 — 7s)/2 and Gn ~ |Ai23| 2 / m f JJ - (Recall: in the SM Csm = £(z/ M <-> v £ ,Gn — > 
Gj?).) Let us now consider the KARMEN experiment 0. Muon anti-neutrinos are produced 
in muon decay, /i + — > e + v s e v s ^ and electron anti-neutrinos are searched through inverse beta 
decay, v^p — > e + n. Since tr couples only to leptons, the detector basis coincides with the 
weak basis, but the source basis is different. Muon decay is mediated by W or tr exchange. 
In the weak basis, the W diagram produces uj^ , while the tr diagram produces . The 
strongest neutrino sector independent bound on G^ is obtained from tests of universality in 
\x and P decays |fT2| , [TT]| . From the lower bound flcifl , J2i=i \V U i\ 2 > 0.995, we get 



G 



2 



sin 2 ^~ (j^) <5x 10- 3 . (18) 

We like to show how the recent 90% CL bound from KARMEN 

P eil < 3.1 x 1(T 3 , (19) 

can be used to set stronger bounds on G^ in part of the neutrino sector parameter space. 
We study two limiting cases. For massless neutrinos, from Eq.([L4]) we get the bound 

^) 2 <3.1xl(T 3 , (20) 

which is stronger than the bound (|T8|). In the large Am 2 limit, from Eq.(|T5|) we get the 
combine bound 

Gn\ 2 _ f G 



#) + 2 (|^)^ + 2^^3.1xl0- 3 . (21) 



In the second example we consider the minimal Left Right Symmetric (LRS) model |14 
In this model there is a Higgs triplet, A^, with the Yukawa couplings to leptons (in the 
weak basis), C = fijLfCiT 2 A L P L Lj, where L { are the lepton doublets and C is the charge 
conjugation matrix. A^ exchange leads to the effective four fermion interaction in Eq.([l7|) 
but with Gn ~ |/n/22|/'^ 2 x +- We again consider the KARMEN experiment. Since A^ 
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couples only to leptons, the detector basis coincides with the weak basis, but the source 
basis is different. Muon decay is mediated by W or A^ exchange. In the weak basis, the W 
diagram produces uYf , while the A \ diagram produces . The strongest neutrino sector 
independent bound on G^ is obtained from tests of universality in /i and (3 decays and is 
given in Eq. (|I8D . Therefore, the bounds ( |20|) and (|2T| ) also hold for the effective interaction 
arising from A^ exchange in the minimal LRS model. 

In the third example we consider models with light leptoquarks (LQ) jnj. There are 
several types of LQ that can lead to sufficiently large new neutrino interaction. We con- 
centrate on the {I$)y = (0)i/3 scalar LQ, S, which couples to fermions (in the weak basis), 
C = XijQ^^PhLiS ', where Qi (Lj) are the quark (lepton) doublets. S exchange leads to the 
effective four fermion interaction (in the weak basis) 

\u^P L d)(u^^P L ^) + (d^P L u)(i?^^P L T)} + h.c. (22) 



£ = G 



N 



with Gn ~ |^2iA3i|/m|g. We assume A32 <C A31. Therefore, LQ interactions involving 
strange quarks are negligible. Let us consider the CHORUS, NOMAD and E803 experiments 
I6| . Muon neutrinos are produced in pion and kaon decays, 7r — > fiu and K — > fiu, and tau 



neutrinos are searched through vn — > rp. The new interaction ( p2|) contributes the same to 
pion decay and to the detection process. Therefore, the neutrino produced in pion decay is 
z/J. This illustrates the above mentioned result: had the neutrino beam been produced only 
from pion decay, we could not probe LQ exchange in those experiments. However, the new 
interactions (|22|) do not contribute to kaon decay and the neutrino produced in K — > \w is 



a weak eigenstate muon neutrino, vf[ ' . The strongest neutrino sector independent bound on 
Gat is obtained from the bound on BR(r — ► 7r yu) 

G N \ 2 _ BR(r->7rV) 

.G F J ~ 2BR(r -+tt-u)' 1 ' 

Using |g BR(r -> 7r°/i) < 4.4 x and BR(r -> iz~v) w 0.117 we get 

sm 2 ^~ <2xl(T 4 . (24) 



We like to show how the expected sensitivity of CHORUS and NOMAD, P^ p ~ 10" 4 and 
E803, P^ p ~ 10 -5 |TB| can be used to probe LQ exchange in part of the neutrino sector 
parameter space. Since we concentrate on LQ that can be tested only by neutrinos from 
kaon decay, we have to use the relative appearance probability of neutrinos from kaon decay, 
where we expect P^ P (K — >■ \w) = few x P^ p . We learn that LQ exchange can be tested, 
probably at CHORUS and NOMAD, and certainly by E803. For example, assuming that 
the bound P^ P (K — > \w) ^ 5 x 10~ 5 will be achieved by E803. For massless neutrinos, 
Eq.(|l4D gives 



2 



^5xl0- 5 , (25) 



which is stronger than the bound (24). In the large Am limit, Eq.([U|) gives 



°-) 2 + 2 (^-) Omw + 26 2 mW < 5 x 10- 5 . (26) 



Gp ) \Gp 
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V. SUMMARY 



There are two important ways in which physics beyond the SM can affect the neutrino 
sector: It may give non- vanishing neutrino masses, and it may modify neutrino interactions. 
In this paper we showed how neutrino oscillation experiments probe both effects, and how 
they can be distinguished in some cases. A distance-independent signal can arise from both 
effects. However, an oscillation pattern can arise only when neutrinos are massive. We study 
the condition on the relative strength of the non-standard neutrino interactions in order that 
it can be tested (see Eq. ([16])). Thus, current experiments aimed to reach a sensitivity of the 
order P^ xp « 10~ 4 JTIJ can typically probe new neutrino interactions arising from physics at 
the 1 TeV scale. There are several well motivated NP scenarios that introduce non-diagonal 
couplings that can be probed in this way. Higgs triplet exchange in Left-Right Symmetric 
models |14| , Light leptoquarks exchange in various models [jn| and super-particles exchange 



in Supersymmetric models without R-parity pH|JrT[ ]. 

Our results are of particular interest in light of the growing evidences for physics be- 
yond the SM in the neutrino sector, in particular, the solar neutrino problem p]Jl7[|, the 
atmospheric neutrino deficit |18| and the recent LSND result [ I9| . Those results cannot be 



simultaneously accounted for in a simple three generation model, but only in models with 
more parameters. Usually, it is suggested that those results are hints to models with an 
extended neutrino sector f20| . Alternatively, it might be that they signal NP in neutrino 



interactions. Such NP can be significant for the MSW solution of the solar neutrino problem 
fl2| , |2T|| , for atmospheric neutrinos [fj^] and, as we have discussed, for accelerator experiments. 
A comprehensive analysis of all these experiments, including possible NP, is needed. 
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